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We explore the use of exact diagonalization methods for solving the self consistent equations of 
the cellular dynamical mean field theory (CDMFT) for the one dimensional regular and extended 
Hubbard models. We investigate the nature of the Mott transition and convergence of the method 
as a function of cluster size as well as the optimal allocation of computational resources between 
bath and "cluster-impurity" sites, with a view to develop a renormalization group method in higher 
dimensions. We assess the performance of the method by comparing results for the Green's functions 
in both the spin density wave (SDW) and charge density wave (CDW) phases with accurate density 
matrix renormalization group (DMRG) calculations. 
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I. INTRODUCTION 

The dynamical mean field theory (DMFT)S has been 
successful in accessing certain aspects of the non- 
perturbative phenomena in strongly correlated electron 
systems. For example, it has given new insights and pre- 
dicted various qualitative trends in the redistribution of 
spectral weight as temperature, pressure and doping are 
varied in optical and photomission experiments, near the 
Mott transition. 

Single site DMFT extracts physical quantities such as 
the self energy from a local impurity model with a self 
consistent bath (also refered to as LISA for local impurity 
self consistent approximation) and by construction this 
approach misses the effects of short range correlations, 
such as the ^-dependence of the self energy, which are 
bound to become important at low temperatures. It also 
does not allow the treatment of phases such as d-wave 
superconductivity, and hence there have beett^i&nificant 
efforts to extend the DMFT methodologyBMfflLlfi. 

Many of the elements of the DMFT method can be 
traced n back to the coherent potential approximation 
(CPA)Ij which has been very successful in disordered sys- 
tems. In this context, the search for extensions of the 
CPA was fraught with many difficulties, and extensions 
of DMFT should be scrutinized carefully in this light. 
The DMFT equations, while having a similar spirit to 
those of CPA for determining the effective medium, are 
considerably more complex because they involve the solu- 
tion of a quantum impurity model, which plays the role of 
the effecive Hamiltonian for the local degrees of freedom 
treated in the DMFT. To analyse these equations new 
concepts and numerous techniques have been developed 
over the last decade. Their incorporation into cluster 
methods is promising but requires substantial new work. 
This paper is a contribution in this direction. In this 
paper we focus on one cluster extension of DMFT, -.the 
Cellular Dynamical Mean Field Theory (CDMFT)tM 
The aim of this paper is to investigate how the exact 
diagonalization aproach which was-so-successful in the 
context of the single site DMFTOiia-Ej, can be used in 



the CDMFT context. 

We solve the self consistent cluster equations for the 
one dimensional Hubbard and extended Hubbard mod- 
els using exact diagonalization of the effective underly- 
ing cluster impurity. We choose these models because 
of their generic nature, as well as the fact that DMFT, 
which becomes exact in the limit of infinite coordination 
numbeio, faces the worst case scenario in one dimen- 
sion. Furthermore in one dimension fairly reliable com- 
putations of static as well as dynamic quantities can be 
made using, t j ho . Density Matrix Renormalization Group 
(DMRG^'tatZl approach to carry out a comparative 
study. 

The exact diagonalization approach, within the single 
site DMFT context, lead to the development of a tech- 
nique inspired by renormalization group ideas which re- 
sulted in the first quantitative ftjidy of the critical prop- 
erties near the Mott transitionta. With a view towards 
the future, we discuss the possibility of combining renor- 
malization group ideas and CDMFT to develop a new 
numerical method in the spirit of the DMRG. 

The Extended Hubbard Model (EHM) is defined by 
the following Hamiltonian: 



H = -*£ ( c W^ + h - c -) + U J2 



UjUj + i — \l 



(i) 



The first term corresponding to hopping between nearest 
neighbor sites and the second term to the onsite Coulomb 
repulsion provide the competition between itineracy and 
localization in the regular Hubbard model. The third 
term represents Coulomb repulsion between electrons oc- 
cupying nearest neighbor sites. The Hamiltonian as writ- 
ten above with /i = U/2 + 2V guarantees an insulating 
ground state with a filling of one electron per site. 

In the following section we introduce and discuss the 
CDMFT self consistency equations and show how they 
can be generalized to treat non-local Coulomb interac- 
tions even in the presence of broken symmetries. In the 
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next and subsequent sections we present results for the 
regular Hubbard and the extended Hubbard model re- 
spectively. This is followed by the conclusions and an 
outlook for further work. 



II. THE CDMFT METHOD 

The CDMFt! method is a straightforward general- 
ization of the single site DMFT, in which local de- 
grees of freedom within a cluster are treated exactly and 
those outside the cluster are replaced by a bath of non- 
interacting electrons determined self consistently. This 
approach, formulated with realistic studies in mind, deals 
in principle with overlapping cells and with clusters which 
are not necessarily defined in real space and could even 
be defined by non-orthogonal orbitals. In this paper, we 
ignore these complications and work with real space clus- 
ters using an orthogonal basis. It is worth pointing out, 
that if we apply the approach to non-interacting disor- 
dered alloys, the self consistency condition of CDMFT 
becomes identical to that iised in the molecular CPA ap- 
proximation of DucastelletJ. This is in the same way as 
the self consistency condition in single site DMFT be- 
comes identical to the single site CPA and the two imptu 
rity approximation proposed_by Ingersent and Schillerc2l 
and by Georges and Kotliaiu, becomes identical to the 
pair CPA. However, unlike the single site case, where the 
transition from real space to momentum space is unam- 
biguous, one has to be careful in interpreting the results 
of the CDMFT equations in fc-space. This is because 
the goal of CDMFT is to obtain the best possible esti- 
mates of local quantities that live within a cluster. For 
this purpose, it introduces a cluster self energy. Long 
distance properties such as the ones contained in the lat- 
tice Green's function and the lattice self energy are then 
inferred from the cluster Green's function or the cluster 
self energy while maintaining causalitytl. 



A. The Cavity Construction 

The most economical way to arrive at the concepts 
of a dynamical mean field theory is via the cavity con- 
struction, which stresses the point made above, that the 
focus of the method is in extracting local quantities. The 
original lattice is divided into equal clusters of size N c . 
Integrating out the degrees of freedom external to a cho- 
sen cluster (labelled 0), we can formally write an effective 
action which would allow the computation of local quan- 
tities as 



ieff 



(2) 



Here j labels individual clusters and [i is an index la- 
belling sites within each cluster. We can split the original 



action into three parts, 

S = S (0) + So + AS, 



(3) 



where S^ ' includes terms outside the chosen cluster (the 
full action with the cluster replaced by a cavity), So com- 
prises of terms solely within the cluster and finally, AS 
includes those terms that couple the cluster with its en- 
vironment. Explicitly, in the case of the EHM, 

ft 3 
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(4) 
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(5) 



Here E 00 includes the hopping matrix as well as the 
chemical potential within the zeroeth cluster and (jfj,, Of) 
denotes all intercluster nearest neighbors. It is to be 
noted that, while the onsite Coulomb interaction U con- 
tributes only to So, a nonlocal interaction such as the 
nearest neighbor Coulomb repulsion V contributes to 
both So as well as AS. This cavity construction is so far 
merely a relabelling of terms in the original action and is 
exact. However approximations need to be made to ac- 
tually access the local properties within the cluster. We 
approximate the effective action of the cluster by keeping 
only the renormalization of the quadratic terms obtained 
after integrating out the degrees of freedom of the sur- 
rounding environment^. Notice that this approximation 
violates the translational symmetry of the original lattice. 
But this is not a problem, since the spirit of the approach 
is to estimate local quantities. The Gaussian approxima- 
tion for the effects of the environment seen by the cluster 
becomes exact in the limit of large coordination^ and 
therefore, our test case in one dimension constitutes the 
worst case scenario. The effective action for the cluster 
can then be approximated as, 



Z 0pa' C 0po 



f dT C lpa \Go *1 c 0up 



Here, the time dependent Weiss field Qq 1 is now a matrix 
in the cluster variables and is a functional of the EHM 
Green's function with the cluster replaced by a cavity. 
Suppressing spin indices, it can be written on the Mat- 
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subara axis as 
Qq 1 (ioj n ) 



itonduv - E au - V 5u V 



fay 



/"• E/ ( n jp')(0) 
(jP'fiu) 

EE ^g<?^(^)^. 

{infill) (jv'fiv) 

Relating the cavity Green's function, G^°\ to the full 
Green's function of the EHM, the following self consis- 
tency condition to determine the Weiss field can be ob- 
tained: 



(7) 



where the local Green's function on the cluster in one 
dimension is defined as (with an obvious generalization 
to higher dimensions) 



Gi oc (iuj n ) 



■k/N c 



dk 



-tt/n c iuj n + (i- E c (iw„) - i(k) 2n/N c 



(8) 



S c is the cluster self energy as obtained from S e g, t(k) 
is the Fourier transform of the cluster hopping matrix 
of the EHM and k is a vector in the reduced Brillouin 
zone of the superlattice ■£-]• Explicitly for a one 

dimensional model with only nearest neighbor hopping, 



iT{k) = -t[<5 M -„ ±1 + ■ 



-ikN c 



e 0-n c + I i-v±i\ 



The functional approach is useful because it clarifies 
the level of approximation used in treating the broken 
symmetry state. In the case that we consider in this 
paper, all the Hartree graphs paa e ad ded to the functional 
(see Eq. 0). The E-DMFTfflEika includes additional 
graphs to the Hartree terms, but the quality of the results 
presented below show that they would make a very small 
contribution. 



C. An Exact Diagonalization Algorithm 

For the purpose of practical implementation of the al- 
gorithm it is convenient to consider the cluster effective 
action to arise from a generalized impurity Hamiltonian, 

H imp = ^2 E ^ C t<7 C ^ + u E n MT n Ml 

pv ka 

+ ]T {VlAa^ + Vk h :A„a ka ) (10) 



Here, t ka represents the dispersion of an auxiliary non- 
interacting bath and V ka „ is the hybridization matrix 
between the bath and the impurity cluster. The second 
and third terms include all the intra-cluster interaction 
terms. In E we lump together the cluster hopping ma- 
trix and chemical potential, as well as the Hartree terms 
(9) arising from non-local interactions between clusters. For 
the EHM 



B. A functional interpretation 

A second approach to derive the above approximation 
is via the construction of a functional of the quantities 
of interest, say the Green's function restricted to a local 
cluster and its supercell translations, such that the sta- 
tionary point of the functional yields those quantitiestiJ. 
In principle such a functional can be constructed order 
by order in a perturbative expansion in the hopping or 
the interactions while constraining the Green's function 
to a fixed value. In practice, sensible approximations to 
the exact functional are contructed by starting with the 
full Baym-Kadanoff functional $[G] and restricting its 
argument from the full Green's function to the Green's 
function defined only inside the cluster and its periodic 
supercell repetitions. This constraint is realized by the 
cluster selLenexgy which plays the role of the Lagrange 
multiplieiElEj'Eij. Differentiation of this functional gives 
the CDMFT equations. Viewed as an approximation to 
the Baym-Kadanoff functional for the full-lattice, this 
procedure may seem a bad approximationcJ. However, 
this construction of the functional should be viewed as 
an approximation to the exact effective action for local 
quantities and CDMFT uses it precisely for such a pur- 
pose. 
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where 



l 3P 



is computed by using the translational in- 



variance of the system at the level of the superlattice. 

In terms of the impurity model the parametrization of 
the Weiss field function is readily obtained as, 



Qq 1 {iuj n ) = iu n - E - A(iu n ); 
where the hybridization function is given by 



{ioj n )[e ka ,Vt S } 



\ * k{T,fj, kcr,v 



k 



£ka 



(12) 



(13) 



Thus the algorithm involves determining the bath pa- 
rameters tka and V k h a self consistently subject to the 
condition given by Eq. |?]. Our implementation re- 
lies on the solution of the impurity model using exact 
diagonalizationO. Starting with an inital guess for the 
bath parameters we obtain the cluster self energy E c , 
and using the self consistency condition determine a new 
Weiss field, (y [ j" 1 (z<x> n ) nGW . In turn, this defines a new 
hybridization function A(ioj n ) ncw . To close the self con- 
sistency loop we project onto a finite subspace of bath 
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size Nb- This projection is carried out using a conjugate 
gradient minimization of the following distance function, 



D 



1 



(n max + 1) N b N c 



E IN 



(14) 



where n max is the number of grid points on the Matsub- 
ara axis. Although we study the CDMFT equations only 
at zero temperature, the self consistency equations are 
solved on the Matsubara axis. 



III. RESULTS FOR THE HUBBARD MODEL 

(V = 0) 



solution exists through all finite values of U, in agree- 
ment with the well known physics of the one dimensional 
Hubbard model. We measure the half spectral gap as the 
point where the strength of the lowest frequency pole falls 
to 9% of its peak height. This percentage is arrived at 
by requiring that the gap at strong coupling be fixed at 
the exact value. We also find, as exnected from the cav- 
ity construction and noted recentljO, that in the case of 
CDMFT, the bath only couples to sites on the boundary 
of the cluster, whereas for other cluster methods such as 
the dynamical cluster approximation (DCA) all the sites 
in the cluster are equivalently coupled to the bath. The 
self consistent procedure is robust and generates these 
two different kinds of solutions regardless of the nature 
of the starting guess. 



In order to highlight the differences between the sin- 
gle impurity (LISA) and cluster dynamical mean field 
schemes we present results obtained for the local spec- 
tral gap as a function of the onsite Coulomb repulsion 
U. To keep the computational cost similar, we initially 
fix the total number of sites in the cluster and bath, 
N s = N c + Nb = 6, across schemes. The exact result 
for the gap is known from Bethe ansatz and is given byt£l 



A(U) 
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(15) 
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FIG. 1: Spectral gap as a function of U/t for the half filled 
Hubbard model. 

Fig. [l] shows that there is a characteristically different 
behavior between the single impurity (N c — 1) and clus- 
ter schemes (N c = 2). In the single impurity case, for 
U > 7.5t the gap follows the exact result, whereas for 
U < 7t it reduces to values much smaller than the exact 
gap and approaches zero. For the range 7t < U < 7.5t we 
observe a coexistence of the gapless and gapful phases. 
This transition from an insulating to a metallic phase 
is an artifact of the mean field approach which incor- 
porates the physics of higher dimensions where the Mott 
transition is indeed present. On the other hand, CDMFT 
compares excellently to the exact gap and an insulating 




FIG. 2: Imaginary part of the onsite Green's function (left 
column) and real part of the nearest neighbor Green's function 
(right column) on the Matsubara axis for different values of 
U (in units of t). 

To further illustrate the performance of the CDMFT 
method, in Fig. || we compare the onsite as well as near- 
est neighbor Green's functions on the Matsubara axis to 
those we obtain using the DMRG method. The calcula- 
tion of dynamical correlation functions within the DMRG. 
method is carried out by using the finite size .algorithmE£l 
combined with the Lanczos vector methodliSO. The re- 
sults shown in the figures are those for chains with 18 
sites and open boundary conditions. We verify that the 
results are close to those found by even longer chains 
and hence representative of the system in the thermody- 
namic limit. To highlight the role of the bath within the 
CDMFT method, we also overlay results obtained from 
exact diagonalization (ED) of a cluster of size N c = 2, 
but without the bath. As U decreases we expect the 
difference between onsite Green's functions from the iso- 
lated cluster (ED) and the DMRG to increase, since the 
physics becomes less local. The contribution of the bath 
in CDMFT therefore becomes more significant with de- 
creasing U and shows a systematic improvement with 
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FIG. 3: Deviation of LISA and CDMFT Green's functions from the DMRG results as a function of the bath size Nh 



bath size as will be discussed later. Throughout the range 
of U we consider, the imaginary part of the onsite Green's 
function computed using CDMFT matches remarkably 
well with the DMRG curve. On the other hand a com- 
parison for the real part of the nearest neighbour Green's 
function shows that the self consistent CDMFT solution 
achieves only a marginal improvement over the ED cal- 
culation for small U and almost none for large U. It can 
be clearly seen that the single impurity approach (LISA) 
performs poorly for the onsite Green's function for all 
U, except when U is large. Likewise, the nearest neigh- 
bor Green's function computed using LISA agrees poorly 
with the exact result throughout the entire range in U. 

An important question to address is the issue of con- 
vergence and the nature of the self consistent solution 
as a function of the size of the effective impurity cluster 
and bath allotted in each case. In order to find a satis- 
factory self consistent solution, we find that the size of 
the Hilbert space allotted to the bath has to be at least 
comparable to that alloted to the cluster itself. For in- 
stance, clusters of size N c — 2 require a minimum bath 
of size Nb = 4 for good convergence. Even and odd clus- 
ters show qualitatively different types of solutions, as we 
discuss further below. In Fig. ^ we plot x, a measure of 
the deviation of the LISA and CDMFT solutions for the 
onsite and nearest neighbor Green's function from accu- 
rate results obtained using DMRG, as a function of the 
bath size Nb. This measure x, is defined as the integrated 
absolute difference between the CDMFT (or LISA) and 
DMRG Green's functions. 

Let us first discuss the results for the strong coupling 
case (U = 12t). The onsite Green's function computed 
using LISA compares excellently with DMRG and this 
continues to be so with the CDMFT method. There is 
practically no dependence of \ on the bath size. For the 
offsite Green's function, x shows a systematic reduction 
with increasing cluster size, but remains fairly indepen- 
dent of the bath size. 

For weak coupling, both LISA and CDMFT become 
exact. The toughest case is when U is of the order of 
the bandwidth; so in the following we discuss the results 
for U — 2t. For even clusters, we find that x f° r the 
onsite Green's function shows a systematic decrease with 
increasing bath size once a minimum bath size is reached. 



On the other hand, x f° r the offsite Green's function 
shows no definite trend with increasing bath size. 
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FIG. 4: From left to right we show the local spectral func- 
tions for sites one through three in a cluster of size N c = 3. 
The bath size is fixed at Nb = 6. 

Contrary to the expectations for a mean field solution, 
even cluster sizes (N c = 2, 4) show no explicit local spin 
symmetry breaking. Odd clusters (N c = 3, 5) show simi- 
lar behavior for small bath sizes upto Nb = N c . For larger 
bath sizes, Nb > N c , and large enough U, the behavior 
of the local spectral function is spin dependent as can be 
seen from Fig. || for U = 12t. Odd clusters with broken 
symmetry are clearly inconsistent with cluster periodicty 
and consequently show poorer convergence. Thus, for all 
even clusters and odd clusters with small enough bath 
sizes, the one dimensional character of the problem is 
dominant over the role of the bath and prevents the self 
consistency from showing explicit spin symmetry break- 
ing. On the other hand, odd clusters (we were able to 
test only the case N c = 3) with large U and a sufficiently 
large bath show a symmetry broken solution consistent 
with the mean field approach. 

The effect of the bath for N c = 3 has an even more 
dramatic consequence at smaller U. As opposed to even 
clusters which correctly see only an insulating solution in 
one dimension, we see both a metallic and an insulating 
solution for Nb > 5 and U < 3t. In the region of coex- 
istence, the metallic solution shows a better convergence 
than the insulating one. Thus, the improvement gained 
in going from LISA to N c = 3 is that the region of coexis- 
tence moves to smaller C/'s and the insulating solution is 
present all the way upto U = 0. Small baths, due to the 
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absence of enough degrees of freedom to impose a mean 
field character to the solution, show an insulating state 
throughout the range in U. Beyond N c — 3 for large 
enough odd clusters we expect to see only an insulating 
solution for all U > 0. 



IV. RESULTS FOR THE EHM (V 0) 




FIG. 5: Imaginary part of the onsite Green's function (left 
column) and real part of the nearest neighbor Green's function 
(right column) on the Matsubara axis for different values of 
V in the SDW phase. 
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FIG. 7: Imaginary part of the onsite Green's function in the 
CDW phase. 
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FIG. 8: Real part of the nearest neighbor Green's function in 
the CDW phase. 
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FIG. 6: Real part of the onsite Green's function in the CDW 
phase. 

In this section we study the performance of cluster 
mean held methods as applied to the Extended Hubbard 
model. This model constitutes a natural test case since 
the non-local nearest neighbor Coulomb repulsion among 
different sites within the cluster can be accounted for ex- 
actly. 

At half filling, for large enough U (strong coupling), the 
system goes through a first order transition from a spin 



density wave (SDW) to a chargerdossity wave (CDW) or- 
dered state at roughly V — U /2c3l2§E2l. We work in this 
regime, by fixing U — 12t and scanning through the near- 
est neighbor Coulomb repulsion V. The CDMFT method 
shows a clear signature of this transition with a cluster 
as small as iV c = 2. In Fig. || we report a comparison 
of the imaginary part of the onsite Green's function and 
the real part of the near neighbor Green's function in the 
SDW phase. As seen in the V — case, the comparison 
with DMRG is better for the former than the latter. We 
note that with increasing V the offsite Green's function 
compares better with the DMRG result. Figs. |]j7| and ^ 
compare the Green's functions within the CDW phase. 
Across the range in V, the agreement with DMRG is so 
remarkable that only one set of curves are discernable 
within the figures. In particular, Fig. |(] shows the real 
part of the onsite Green's function which is zero in the 
SDW phase but acquires a non-zero value in the CDW 
phase. This indicates a breaking of local particle-hole 
symmetry of the system that is now conserved only on 
the average every two sites. Notice how the curves for 
different values of V cross, defining a characteristic scale 
approximately equal to U. The versatility of CDMFT 
in being able to treat the ordered phase stems from the 
fact that it only assumes supercell periodicity and con- 
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sequently allows for complicated ordered states within a 
cell. This is in contrast to the DCA method that assumes 
full periodicity of the lattice, and therefore requires a fine 
discretization in momentum space to adequately capture 
the discontinuities in the Brillouin zone that appear with 
the emergence of short range ordered phases. 



V. CONCLUSIONS AND OUTLOOK 




FIG. 9: Schematic representation of an RG procedure in- 
tended to reach large cluster sizes. 

We have shown that the CDMFT method produces 
remarkably good results when tested in one dimensional 
systems. This is very encouraging, since one would ex- 
pect mean field methods to perform even better as the 
dimensionality increases. We saw that even the smallest 
possible cluster size (N c — 2) along with a self consistent 
bath of modest size allows for an accurate determination 
of onsite correlations, while the offsitc quantities improve 
systematically with increasing cluster size. Further, the 
success of our extension of the CDMFT equations to de- 
scribe the discrete broken symmetry of the CDW order 
clearly demonstrates the power of the method to incor- 
porate short range correlations. 

Following the example of the EHM, with the simplest 
kind of short range correlation, it should be possible to 
treat the physics of more complicated ordered phases as 
well as unit cells of realistic materials within the cellular 
mean field approach. The main obstacle in this program 
is the computational effort that goes into solving the self 
consistent equations for the large clusters required by real 
world applications. This calls for a systematic approach 



that blends together clustep-dypamical mean field and 
renormalization group ideasE3E2l to select the most rele- 
vant degrees of freedom to represent a cluster. 

We suggest that the CDMFT approach can be used 
iteratively in complete analogy with the density matrix 
renormalization group approach. It is useful to compare 
the similarities and differences between these two meth- 
ods. In both of them, one improves on the exact diago- 
nalization of a small finite size system by embedding the 
small system in a larger one; an embedding which can be 
described by a reduced density matrix. Within DMRG 
the density matrix is determined by an exact diagonaliza- 
tion of the larger system containing the subsystem of in- 
terest. In CDMFT, one can approximate the density ma- 
trix by modelling the enviroment with a Gaussian Weiss 
field. 

The CDMFT assumption of a self consistent gaussian 
bath as the environment, that one can integrate out ex- 
actly in order to define a reduced density matrix for the 
cluster, becomes more and more accurate as the dimen- 
sionality increases. This choice of bath is optimum in 
a dynamical mean field sense (optimized for the compu- 
tation of the local one particle Green's function of the 
cluster) and sidesteps the conventional DMRG proce- 
dure for building the reduced density matrix out .-of a 
few target states. Further, the DMRG prescriptiorJia of 
selecting the states for which the reduced density ma- 
trix has the largest eigenvalues allows for a truncation 
of the cluster Hilbert space while retaining the relevant 
physical information. In Fig. ^|we indicate schematically 
how such a method would proceed by doing CDMFT on 
small clusters, truncating their Hilbert space using the 
DMRG prescription and finally using them as building 
blocks of even larger clusters. This procedure can be re- 
peated until either a desired cluster size or convergence 
in a certain observable of interest is reached. We believe 
this scheme should open new vistas for numerical renor- 
malization group calculations of realistic systems in two 
and three dimensions. 
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APPENDIX A: DETERMINATION OF THE 
CLUSTER SELF ENERGY 

The first step in the CDMFT iterative scheme is to 
compute the cluster self energy S c for an initial guess of 
the bath parameters. This can be done using Eq. ^ by 
subtracting the inverses of the numerically determined 
cluster Green's function {Gi mp (iu) n )) from the exactly 
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known Weiss field. Although this is accurate enough for 
most purposes, when the results are small, it can intro- 
duce numerical errors. To deal with this, an alternative 
procedure was proposedJn the context of the single im- 
purity Anderson modelEil. The idea is to isolate contri- 
butions to the self energy coming purely from the hy- 
bridization and the interactions. The terms arising from 
the interactions can be written as ratios of correlation 
functions of composite operators leading to a more sta- 
ble numerical procedure. We generalize this procedure 
to the cluster Anderson impurity model to write the self 
energy in the form 



and 



t c (iu> n ) = £u(iuj n ) + t, v (iu n ) 



where 



imp 



imp 



(Al) 

" (iw n ) (A2) 
(tw„) (A3) 



* a 



(A4) 



(A5) 



Here < J 4,B> denotes the Green's function of operators 
A and B. This method of computing the self energy is 
particularly robust for small values of U/t, when the self 
energy is relatively small in magnitude, as compared to 
the direct approach. 
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